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Elizabeth Gasparim and Irena Swanson * 

■ We present an algorithm for computing instanton numbers of curve singu- 
I larities. A comparison is made between these and some other invariants of 

psj ' curve singularities. The algorithm is implemented in Macaulay2, and can be 

. downloaded from ^ttp: / /www.math.nmsu.edu/~iswanson/instanton.m2| . 

<^ : 

^ , 

Q>^ ■ We start with any polynomial y) defining a plane curve with singularity at G C 

^ I Let 7r: —>■ denote the blow-up of at the origin and let j be a positive integer. 

The data {j,p) determines a holomorphic bundle E{j,p) on with splitting type j and 
^ ! extension class p. We algorithmically compute numerical invariants of the bundle E{j,p) 
r-i ! and use them as invariants of the curve. 

. Tables giving examples of these new invariants of the curve together with some classical 

S ■ invariants are given in Section 1. A specially interesting example appears on table IV, 
where we give two inequivalent singularities that have all the same classical invariants, but 
^ ■ are distinguished by instanton numbers. 

^ ■ The holomorphic bundle E{j^p) and its numerical invariants have interpretation in 

O ' mathematical physics as instantons on and numerical invariants of the instantons. Here 
O . we mention briefly some properties of these invariants. More details are given in [G3]. 
^ ! Instantons are well known to have a topological invariant called the charge, which in the 

r-| ■ case of compact surfaces corresponds to a second Chern number of a vector bundle. For 
. the case of bundles on a blown-up surface, the local Chern number around the exceptional 

S ■ divisor decomposes as a sum of two numerical invariants (denoted w and h), which are 

■ local analytic invariants of the bundle. In [BG] it is shown that the pair of invariants 
^ ' (w, /i) is strictly finer than the Chern numbers in the following sense. The pair (w, h) 

' provides the coarsest stratification of moduli of instantons on the blown-up plane for which 
the strata are Hausdorff. In contrast, the stratification by topological charge does not 
provide Hausdorff strata [BG, Theorem 4.1]. Hence, the pair (w, h) gives strong numerical 
invariants, detecting more than topological information. 

The idea of using these invariants for curves is natural given that, as a first step 
toward resolving the singularity of the plane curve defined by p, one blows up the plane, 
thus arriving at C , the base space of the bundles we construct. 



* The first author acknowledges support from NSF NMSU ADVANCE grant NSF-0123690 and the 
second author acknowledges support from NSF grant DMS-9970566. 
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Firstly, let us explain how to pass from curves on to vector bundles on C^. In 
what follows we fix, once and for all, the following coordinate charts. Let {x,y) be the 
coordinates on C . Write C =U UV where 

U = {iz,u)}^C^^{i^,v)} = V 

with (^, v) = {z~^, zu) inU r\V . Note that in these coordinates, the blow-up map tt: 
gives X = u, y = zu. We denote by £ the exceptional divisor. 
Once these charts are fixed, there is a canonical choice of transition matrix for bundles 
on C^. If i? is a holomorphic rank 2 bundle on with ci{E) = then there exists an integer 
i, called the splitting type of such that E\i ^ 0{j) © 0{-j). By [G2, Theorem 2.1] E 
has a canonical transition matrix of the form 

' q 



(1) 



from U to V, where 



2j-2 j-1 



i=i 

is a polynomial in z, z~^ and u. 

In this paper, we start with a polynomial p e C[a;, y] and for each positive integer j 
we associate to the pair a holomorphic bundle E{j,p) on obtained as follows. 

Definition 0.1: Given a polynomial p{x,y) defined over C[x,y] we first consider the 
polynomial p{u, zu) = YliiYliiPii^^'^^ obtained from p by making x = u and y = zu. If 
furthermore a positive integer j is chosen, we define p = ^iPuz^u'^ by the rule 

p. =h^i if l<i<2j -2and0<l<j -1; 
* [ otherwise. 

We then set E :— E{j,p) to be the bundle given in canonical coordinates (as above) by the 
transition matrix 

We have therefore associated to each pair (j, p) formed by an integer and a polynomial 
p{x,y) a unique vector bundle E(j,p) over C^. We now define instanton invariants w and 
h for these vector bundles. 

Definition 0.2: Given a bundle E over we define the sheaf Q by the exact sequence 

O^TT.E^ (tt.E)^^ ^ Q ^ 
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and set 



w{E): = 1{Q), h{E): = 1{R^tt^E). 



In instanton terminology, w is called the width of the instanton and h is called the 
height of the instanton. The charge of the instanton is c = w+h. These are clearly analytic 
invariants of the bundle. To use them as invariants of the singularity, and to perform the 
calculations, we choose charts for the bundle and a representative p of the singularity. 
Such choice works particularly well in the case of quasi-homogeneous singularities, where 
preferred representatives are known to exist, cf. Saito [Sa] or Arnold [Ar, page 95]. 

If a different representative p' of the germ of the singularity is chosen, then by classical 
theory it is known that there exists a holomorphic change of coordinates taking p to p'. 
The same change of coordinates has to be made to take the bundle E{j^p) to a bundle 
E'{j,p'). E and E' are isomorphic bundles and therefore have the same invariants. In this 
sense the numbers can be considered as analytic invariants for germs of curve singularities. 

In this paper we prove that the computation of h and w can be automated. We 
implemented our algorithm in the symbolic computer algebra package Macaulay2 (due to 
Grayson-Stillman [GS]). Having implemented algorithms to compute both invariants, we 
obtained a large amount of examples, which lead us to a simple formula for h. 

Theorem 3.3: Let m denote the largest power of u dividing p. If E is the bundle defined 
by data {j,p), then h{E) = g) - {^-'^). 

The algorithm we present here has as input data the polynomial p{x, y) and the integer 
J, and as output the values of the instanton height h{E{j,p)) and width w{E{j,p)). Our 
implementation on Macaulay2 is for p with rational coefficients. 

1. Some tables 

Before giving the algorithm, we illustrate some instanton numbers by tabulating them 
together with some classical invariants. Tables I and II give examples of instanton numbers 
of monomials. 



TABLE I 


splitting type 2 


polynomial 


w 


h 


charge 


X 


1 


1 


2 


y 


1 


1 


2 


xy 


2 


1 


3 
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TABLE II 


splitting type 3 


polynomial 


w 


h 


charge 


X 


1 


2 


3 


y 


1 


2 


3 


x^ 


3 


3 


6 


xy 


2 


3 


5 


y' 


3 


3 


6 


x^y 


4 


3 


7 


2 

xy 


4 


3 


7 


x'^y'^ 


5 


3 


8 



Table III illustrates that for a fixed polynomial, we may get different invariants as we vary 
the splitting type. 



TABLE III 


splitting type 7 


splitting type 8 


polynomial 


w 


h 


charge 


w 


h 


charge 


2 2 

x^ — y^ 


2 


11 


13 


3 


13 


16 


— y^ 


3 


11 


U 


3 


13 


16 


x'^ — y^ 


3 


11 


U 


3 


13 


16 


X — y 


4 


15 


19 


4 


18 


22 


3 4 

X — y 


6 


15 


21 


6 


18 


24 


X — y 


6 


15 


21 


6 


18 


24 



Our interest in instanton numbers was partially fueled by the fact that in some cases, 
instanton numbers give finer information than the classical invariants. We considered the 
invariants: multiplicity m, Sp = dim (O/O), = Milnor number, and r = Tjurina number. 
Note that the Milnor and Tjurina numbers are defined only for isolated singularities, but 
instanton numbers are well defined for non-isolated singularities as well. Here is an example 
where instanton numbers distinguish the singularities, but other invariants do not. 



TABLE IV 


splitting type 3 


polynomial 


m 


Sp 




T 


w 


h 


charge 


x^ — x^y + y^ 


3 


3 


4 


4 


4 


3 


7 


x^ — x'^y'^ + y^ 


3 


3 


4 


4 


5 


3 


8 
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Note that the first polynomial is irreducible, whereas the second is reducible in the local 
ring, cf. [Ha, ex. I 5.14], so they define inequivalent singularities. 

We believe that the fact that the classical invariants m, 5p, /i, and r do not distinguish 
some inequivalent singularities is an evidence that finer invariants are useful. In tables 
V-VII below we give examples of inequivalent singularities which are distinguished by 
instanton numbers, and also by one other classical invariant. 



TABLE V 


splitting type 4 


polynomial 


m 


Sp 


fi 


T 


w 


h 


charge 


— y'^ 


2 


3 


6 


6 


3 


5 


8 


— y^ 


3 


3 


6 


6 


6 


6 


12 



Table VI comes from [Ha, ex. V 3.8] on page 395. However, in the statement of this 
exercise, the first polynomial contains an incorrect exponent: it is written as "x^y — y^" 
but it should be "x^y — y^" . 



TABLE VI 


splitting type 4 


polynomial 


in 


dp 


II 


T 


a: 


h 


charge 


4 5 

X — xy 


4 


9 


17 


17 


10 


6 


16 


x^ — x^y^ - x^y^ - y^ 


4 


9 


17 


15 


8 


6 


14 



Table VII comes from the list of bimodular singularities given by Arnold in [Ar, page 159]. 



TABLE VII 


splitting type 8 


polynomial 


m 


Sp 


fi 


T 


w 


h 


charge 


x^ + x'^y^ + y^ + xy'^ 


3 


9 


15 


16 


6 


18 


24 


x^y + x'^y^ + xy^ + y'^ 


4 


9 


U 


15 


7 


22 


29 


x^ + x'^y^ + y^ 


4 


8 


15 


15 


10 


22 


32 




4 


7 


12 


13 


9 


22 


31 


{x'^ + y^Y + ^y^ 


4 


6 


9 


9 


6 


22 


28 



Table VIII below comes from the list of exceptional families of unimodal singularities in 
[Ar, page 95]. 
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TABLE VIII 

-1. ^ 1. J—/ J— J A 1 \ J. J. J. 


splitting type 7 


pUiy 1 LU 1 1 ibthl 




11 


charge 







1 K 




_1_ 7/8 _1_ XI 


6 


15 


21 




7 


18 


25 


+ xy^ + 




15 


21 


a^^y + + a^Z/^ 


7 


18 


25 




7 


18 


25 


4 1 4 1 6 




18 


25 



2. Computing the instanton width 

In this section we present an algorithm that takes as input the pair and computes 
the instanton width w = 1{Q), where Q is the skyscraper sheaf as given in Definition 0.2. 
Since Q is supported at zero, 1{Q) is the dimension of the stalk at zero. Let M be the 
completion of the stalk 7r*-E at zero, that is, M:= (7r*-E)o. Then the length of Q equals 
the dimension of the cokernel of the canonical map from M to its double dual. If we can 
compute M, we can also compute Q via the following lemma: 

Lemma 2.1: Let R be a commutative Noetherian ring and A an n x m matrix with 
entries in R. Let M he the R-module such that R^ — >R^ — > M — > zs an exact sequence. 
Let N be the kernel of the transpose of A. Then N is a submodule of R^, say generated by 
t elements. Let B be the n x t matrix whose columns are the given generators of N. Let 
C be the matrix such that R^ — >R^ N ^ is exact. Let _^ = Homn{_ ,R). Then 
(i) M^^ is isomorphic to the kernel of the transpose of C , 

(a) the image of the canonical map M — >• M^^ is isomorphic to the module generated by 

the columns of the transpose of B, 
(Hi) M'^'^ /M is isomorphic to {kernel C'^) / {image -B^). 

Proof: By left-exactness of the Horn functor, = HomR{M, R) is the kernel of the 
map defined by the transpose of A. Thus = M^. Similarly, by the definition of C, M^^ 
is the kernel of the transpose of C. Let /i, . . . , /t generate C R^. The natural map 
M — > M^^ takes m to (/j^, where for each / G M^, ipm{f) = /("^)- Thus, in coordinates, 
the image of m G M in M^^ C R^ is the vector (/i(m), . . . , ft{m)), therefore the image of 
M in M^^ is generated by the columns of the transpose of B. D 

Once Q is computed as in the previous lemma, its length can be computed as well. 
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Thus it remains to compute the C[[a;, j/]] -module M. By the Theorem on Formal Functions, 

M^\mlH\£^,E\er^) 

as C[[a;, |/]]-modules, where is the nth formal neighborhood of £. We use the following 
lemma: 

Lemma 2.2: [Gl, Lemma 2.2] Set M = H^{i2j-2, E\e^._^), and let p: M ^ M^^ be the 
natural inclusion of M into its bidual. Then 1{Q) = dim cokerp. 

This lemma greatly simplifies the calculations. Using this lemma, for the purpose of 
finding 1{Q), we may assume M = M. Under this identification, and with our choice of 
transition matrix T in (3), it then follows that M is generated by all holomorphic vectors 
V over C[[2;, tt]] for which Tv is a holomorphic vector over C[[2;~^, zu]]. Thus, we need to 
find all vectors v = (a, b), given by a = i b = ; buz^u^ , where all i and I are 

non-negative, and au^bu are (unknowns) in the field, for which Tv is holomorphic over 
C[[2~^, This restriction yields relations among the unknowns aujbu. To get all such 
relations, at the beginning we treat an and bu as variables. 

Lemma 2.3: Whenever I > i + j, then bu = 0. Whenever I > i, then an = 0. 

Proof: The second coordinate of Tv is z~^b = ^ - ^ biiz^~^u^ = ^ biiz^~^~'^{zuy . In order 
for this to be holomorphic in z~^ and zu, necessarily the coefficients bu with I > i+j must 
vanish. 

The first entry of Tv is 

z^ auz^u' +pY^ bi^iy^z^^ 

2j-2 j-1 

i,l io=l lo—io—j+'i- hjh 

Let m be the minimum ?i-dcgrec occurring in p, i.e., the minimum zq such that p^,,;,, ^ 0. 
Then for each i < the coefficient of in the first entry of Tv is Yli au^zuY z''~^^K In 
order for this to be holomorphic in zu and z~^ , necessarily for all / > i — j, au ~ 0. Now 
consider the case i > m. The coefficient of {zuyz^~^^~'^ in the first entry of Tv is 

2j-2 j-1 
io = l io=io— j + 1 

where bi>^i> is treated as zero whenever i' or I' is negative. By the established bounds, 
whenever I > i — io + Iq, then bi-i^^i+j-ii^ = 0. Since p is a subpolynomial of p{u, zu), the 
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only pairs (zq, ^o) to consider are those with Iq < ig. Thus whenever I > i — io + Iq, then 
bi-ig^i^j^i^^ = 0. In particular, if I > i, then all are zero, so that the coefficient 

of {zuyz^~^^~^ in the first entry of Tv is an. But if Tv is to be holomorphic in zu,z~^, 
then necessarily an = 0. D 

Thus for each i, we need only to consider finitely many unknowns an, bu to construct 
the vectors v as in the set-up above, and find relations on these. Using Lemma 2.2, we may 
assume that M = H^{l2j-2iE\(^^._^). This means that we only need relations involving 
an, bu for z < 2_7 — 2. But the relations involving bu with z < 2_7 — 2 arising from z^a + pb 
being holomorphic over C[[tt, zu]] involve variables aiq and bi'i with z' < 2_7 — 2 + degup, 
so that 

Lemma 2.4: Let N be the sum of 2j — 2 plus the u-degree of p. Then M is generated 
by vectors (a, b) with a = ^ auz^u^ , b = ^ buz^u^ , where all i and I are non-negative, 
au,bu e C, and 

(i) au = bu = whenever i > N. 
(a) bu = whenever I > i + j. 

(Hi) au — whenever I > i. D 

The conditions on the au and bu as in the lemma above are not the only ones that 
arise from the condition that Tv be holomorphic in z~^,zu. Finding all the conditions 
amounts to finding the generators of M. 

Definition 2.5: All the remaining relations on the variables au, bu arise from the condition 
that in the first entry ofTv, whenever I > i, then the coefficient of u^z^ must be 0. From 
now on, we refer to these coefficients as the generating relations. 

Note that these coefficients are all linear forms in the ring 'C\_au, bii\. 

As in [Gl], we find these relations successively in the zeroth through A^th neighbor- 
hoods. Below follows our algorithm which computes the generating relations, arising from 
ensuring that the first coordinate of Tv is holomorphic in zu, z~^. We write all algorithms 
in this paper in pseudocode, close to the Macaulay2 code that we implemented. 

Algorithm "getrelations" to get the generating relations on the an, bu 
Input: non- negative integer N, 
fTv = first entry of Tv 

ring R, the polynomial ring in u, z, all an, bn 
Output: ideals nonfree and relations in R 

nonfree = zero ideal in R 
relations = zero ideal in R 
k = 

while (k < N) ( 

tempoly = truncation of fTv to terms of tt-degree at most k 
while (tempoly != 0) ( 
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tempterm = leading term of tempoly 
i = u exponent of tempterm 
1 = z exponent of tempterm 

partp = coefficient of u*z' in tempoly, linear form in am, hu 
tempoly = tempoly - partp * u*z' 

fTv = fTv - partp * u*z^ 

if (1 > i) then ( 

relations = relations + ideal (partp) 

nonfree = nonfree + ideal (leading variable in partp) 

) 

) 

k = k + 1 

) 

return nonfree, relations 

) 

The ideal of all the generating relations obtained in this way is called relations. The 
ideal nonfree contains the leading variable of each relation: in the sense of linear algebra, 
these leading variables are not free. We use the following ordering of the monomials. 

Definition 2.6: We say that an > ai'i' if i > i' or if i = i' and I > V , and similarly 
bii > bi'i' if i > i' or if i = i' and I > V , and furthermore that an > bi'i' for all i, I, i', V . 

Thus by the form of Tv, each relation contains at most one an, and each an appears 
at most once in a generating relation. Observe that if all the coefficients of p{x,y) are in 
a subfield F of C, then all the generating relations have coefficients in F. The algorithm 
"getrelations" above computes the relations among the given an,bn, but some of these 
relations are "fake" in the following sense: 

Example: Let p{x,y) = — y^,j — 3. Then with notations as above, N = 7, p = 
V? —u^z^. For holomorphic a = ^ ai^u^^;^, 6 = ^ bnu^z'', the coefficient of u^z^ in z^a+pb 
is = asG + bQQ — b^Q. However, if we restrict the first index of an to only vary from to 
N, then "getrelations" gives the "fake" relation = ^69 — ^56- 

Thus, computation of the instanton width will have to account for and remove such 
"fake" relations. We do this as follows. Such relations only involve variables an,bn with 
i > 2j — 2. Thus these variables are not allowed to be free variables in the sense of linear 
algebra. Using Lemma 2.2, the remaining free variables do give a generating set of M as 
follows. For each of the free variables, set that variable to 1 and all the others to in 
(a, b). This produces a finite generating set of M as a module over C[[x, y]] — C[[u, zu]] of 
elements whose entries are in C[[w, z]]. As tt is a non-zero divisor, M is isomorphic to u'^M 
for arbitrary integer k. By Lemma 2.4, u^a and u^b are both polynomials in u{= x) and 
zu{— y), so that the generators of M can be written as pairs of polynomials in C[[x, y]]. 

Algorithm "polyconv" to convert C[w, 2;] -polynomials to C[x, y] -polynomials 
Input: a polynomial f in u and z 
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Output: a "truncated" polynomial g(x, y) such that g(u,zu) = f, 
where f is that part of f for which this can be done 

g = zero element of C[u,z] 
while (f != 0) ( 

If = leading term of f 

i := the u exponent of If 

1 := the z exponent of If 

if (1 < i) then 

g = g + y' * x*~' * leading coefficient of If 

f = f - If 

) 

return g 
) 

Algorithm "setvectors" to express generators of u^~^^ M 

as vectors with entries in C[a;, y] 
Input: polynomials u-" a, u-" b, named Apoly, Bpoly, respectively, 

lists changeables, allvars 
Output: A, a presentation matrix for the C[[x, yjj-module M 

Mxy = zero submodule of C[x,y]'^ 
total = #changeables 

k = 

while (k < total) ( 

tapoly = substitute in Apoly the k*'* changeable variable to 1 

tbpoly = substitute in Bpoly the changeable variable to 1 
tapoly = substitute in tapoly all other variables to 
tbpoly = substitute in tbpoly all other variables to 
A = convert tapoly into a polynomial in x and y (use polyconv) 
B = convert tbpoly into a polynomial in x and y (use polyconv) 

Mxy = Mxy + submodule of C[x,y]^ generated by (A,B) 
k = k + 1 

) 

return a presentation of the C[x, |/] -module Mxy 

) 

The output of the last routine is the presentation matrix of a C[x, y]-module, which 
by faithful flatness of C[[x, y]] over C[x, y] is also the presentation matrix of the C[[x, y]]- 
module M. Finally, tying it all together: 

Algorithm to compute instanton width of the instanton with data (j,p) 
Input: a polynomial p in C[2;, and a non-negative integer j 
Output: the width of the instanton with data (j,p) 

p = p(u,zu) truncated to u-degree at most 2j-2 

N = 2j-2 + u-degree of p 

R = C[u,z,ai/, hii'\ , i < N+1, ordered as in Definition 2.6 

a = ^ii 

b = bii 

fTv = z-' a + p b 

compute relations and nonfree variables as in algorithm getrelations 
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changeables = all a^i , hu with i < 2j — 2, which are not in nonfree 

a = a after applying all the relations 

b = b after applying all the relations 

A the presenting matrix of the C [x,y] -module M, 

output of setvectorsCu-'+'^apoly, u-^+^bpoly , changeables, allvars) 
Q = cokernel of the natural map M — > M^^ , as in Lemma 2 . 1 
return (length of Q) 

) 

3. Computing the instanton height 

In this section we compute the instanton height. Recall that the instanton height h 
is the length l{E}'K^E{j^p)). Another use of the Theorem on Formal Functions gives 

(i2V£;)^ = lim H\£n,EUJ. 

We use the following lemma: 

Lemma 3.1: [Gl, Lemma 2.3] 1{R^7T^E) = dim£ H^{£2j-2, E\i^._^). 

This lemma greatly simplifies the calculations. We then proceed to compute the first 
Cech cohomology. 

Remark: We have two options to compute Cech cohomology. Given that our open 
cover of the base space is given by affine sets that are open in the analytic topology as 
well as in the Zariski topology, we have the option to compute either holomorphic Cech 
cohomology (by taking holomorphic cochains) or else algebraic Cech cohomology (by taking 
algebraic, i.e, polynomial cochains). Since £ is compact, so are its formal neighborhoods. 
By Serre's G.A.G.A., holomorphic bundles on a compact variety are algebraic, therefore 
Hlig{£n, E\£^) = Hl^i{£n, E\£^). It follows that the holomorphic and the algebraic methods 
give the same answer. 

We compute the instanton height using holomorphic Cech cohomology. The 1-cocycles 
consist of the vectors (a, b) which are holomorphic functions defined on the intersection 
U nV. Hence a,b E C[[u, z, z~^]]. The coboundaries consist of the vectors of the form 
V + T~^v', where v is holomorphic in z, u (on U) and v' in z~^ .,zu (on V). First of all we 
choose simple representatives for the cocycles: 

Lemma 3.2: Every 1-cocycle has a representative of the form 

i=0 i=i-j+i \ / 

with an e C. In particular, every 1-cochain represented by ^'^q^^ ^^"^^ i,l > Q is a 
coboundary. 
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Proof: Let u be a 1-cocycle and let ~ denote cohomological equivalence. A power series 
representative for a l-cochain has the form a = 7 7 ( ] ^'w*, with au^hu e C. 

The l-cochain si — ( , J is holomorphic in [/, hence is a coboundary. Hence 

i=o 1=0 V »^ / 



i=0 i=— 00 



After a change of coordinates 

00 

^-E E 



i=0 l= — oo 



z-%1 ' " ^ ' 



but given that S2 — f j z^u^ is holomorphic in 



i=0 Z= — 00 

CX) —1 

Tcr ~ Ta — S2 — 

1=0 l= — oo 

and going back to the [/-coordinate chart, 

00 — 1 



z^au+phii \ I 



i=0 l= — oo 



au + z ^pbii \ I i 

• 



But p contains only terms 2;^^ for A; < j — 1, therefore z ^p contains only negative powers of 
z. Renaming the coefficients we may write a — ( / ^^^^ some a^^ G C, and 

i=0 l=-oo ^ ^ 

consequently Tcr = ( Here each term a^^^-^+^it* satisfying j + I < i 

i=0 l=-oo ^ ^ 

is holomorphic in the \^-chart. Subtracting these holomorphic terms we are left with an 
expression for a where the index I varies asz— j-|-l<i<— 1. This in turn forces i < j — 2, 
giving the claimed expression for the 1-cocycle. D 

Thus we only have to consider cocycles of the form (a, 0). Which of the cocycles (a, 0) 
is a coboundary? In other words, which (a, 0) equal to (c, d) + T~^{c', d'), where c, d are 
holomorphic on U and c',d' on V, or even more simply, for what c and d holomorphic on 
U is T(a + c, d) holomorphic on V? The second coordinate of T(a + c, d) is z~^d, and in 
order for that to be holomorphic, du = whenever I > i + j. This is the only restriction 
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on c and d obtained from the second coordinate. From the first coordinate of T{a + c, d) 
we obtain the constraint that z^{a + c) +pd be holomorphic on V, that is, holomorphic on 
coordinates and zu. 



Lemma 3.3: Let E he the bundle defined by data {j,p), and m the smallest exponent of 
u appearing in p. Then 



1{R^7T^E) > 



Proof: Let a — (a, 0) denote a 1-cocycle where a = z^u^ with < i < m — 1 and 
i — j + 1 < I < —1. We claim that a represents a nonzero cohomology class. In fact, for 
cj to be a coboundary there must exist c and d, holomorphic in [/, making the expression 
{a + c) + holomorphic in V . However, a is not holomorphic in V. Moreover, by 
the choice of m, no term in pd cancels z^ a. Consequently, no choice of c and d solves the 
problem of holomorphicity on V . Hence l{R^n^E) is at least the number of independent 
cocycles of the form a — (a, 0), where a — z^u^ with < z < rn — 1 and i — j + 1 <l < —1- 
There are (^) - {^~^) such terms. □ 

Theorem 3.4: Let E be the bundle defined by data {j,p), and m the smallest exponent 
of u appearing in p. Then 

Proof: First note that, by Lemma 3.2, if I > j, then (2;'~-'w*,0) is a coboundary. Using 
the proof of Lemma 3.3, it suffices to show that ii j > I > i > m, then a = (2;'~^tt\ 0) 
is a coboundary. Write p = Ylt=o cirz'^~^u'^ + p' where s e {0, . . . , m} is some integer, 
ttr are constants, ag ^ 0, and p' is a polynomial in u and zu each of whose terms has 
■u— degree at least m + 1. Observe that d = a~^u'^~'^z^~'^^^ is holomorphic on C/, since 
by assumption i > m and I > m — s. Therefore a a' — {z''~^u^, —d), where ~ denotes 
cohomological equivalence. Changing coordinates we have Ta' = {z^u^ —pd, —z~^d). From 
j > I > i > m > s we deduce that z~^d is holomorphic on V. We rewrite the first entry 
Ta' as 

s-l 

z\' -pd = J]a;^a^^'+"-''«* + a7^«*-"^^'-"^+V. 

r=0 

By assumption, I— j-Vs — r<s — r<s<m<%^ therefore the first sum on the right side 
of this expression is holomorphic on V . It follows that Ta' ~ {a~^u^~'^ z^~'^^^p' 

Let z^u^ be an arbitrary term in p' . Then r > m and v < r. If i — m + r > 
I — m + s + V, then the term c = aj^u^~'^z^~'^'^^ ■ u^z'" is holomorphic on F, therefore 
Ta' ~ (a^^M^-^^'-^^+V - c, 0). 
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Removing all such terms c, we may now write Ta 0), where 

p contains only terms in z'"u'^ such that i — m -\- r < I — m + s + v. Consequently a' ~ 
{z~^ a~^u'^~^ z^~'^^^ p^ 0) and as i — m + r > z, by (reverse) induction on i and /, each term 
[u'l'-T^ ■ z^ z~^ is a coboundary. Hence, a' is a sum of coboundaries, and is 
itself a coboundary. D 

Thus when starting with p{x,y) G C[x,y], the computation of the instanton height 
is very fast: once m is determined, then the following routine "iheight" finishes the 
computation: 

Algorithm "iheight" to compute instanton height 
Input: a polynomial p and a non-negative integer j 
Output: returns the length of R^. 

p = p(u,zu) truncated to u-degree at most 2j-2 
m = the largest power of u dividing p 
M = j*(j-l)/2; 

if (j > m+1) then M = M - (j-in)*(j-in-l)/2; 
return M 

) 

We implemented in Macaulay2 these algorithms for computing the instanton widths 
and heights. The computation of height of course only takes a negligible amount of time, 
and the computation of widths takes a few seconds. For example, "iwidth(x"4-x*y"5 ,4) " 
finishes in a Linux workstation in 17.07 seconds, and "iwidth(x~4-x~2*y~3 -x~3*y~5- 
y"8,4)" finishes in 32.02 seconds. 

Acknowledgment. We thank Michael Stillman for help with Macaulay2 code. We also 
thank the referees for suggesting several improvements to the exposition. 
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